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1 Abstract 

We study the inverse problem for the so-called operators with energy 
depending potentials. In particular, we study spectral operators with 
quadratic dependance on the spectral parameter. The corresponding 
hierarchy of integrable equations includes the Kaup-Bousinesq equa- 
C/3 | tion. We formulate the inverse problem as a Riemann-Hilbert prob- 

£^ ■ lem with a Z2 reduction group. The soliton solutions are explicitly 

obtained. 
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1 Introduction 

The last decades witnessed an explosion in the complexity and sophistication 
of mathematical theories for fluids and in particular for water waves. The soli- 
ton theory has been always at the center of these developments, such as from 
its early days the soliton theory has transformed and enhanced enormously 
the mathematical description of nonlinear wave propagation. The simplest 
and best known integrable water-wave equations belong to the Korteweg-de 
Vries family. For some classical and modern aspects of the theory of water 
waves, nonlinear waves and soliton theory we refer to the following mono- 
graphs and the references therein: rjl[3[9l[IOl[ia[IS[]Sl[2a 

There are classes of soliton equations whose associated spectral problems 
are polynomial in the spectral parameter. They are known also as soliton 



1 E-mail: Rossen. Ivanov@dit.ie, 
2 E-mail: Tony.Lyons@mydit.ie, 



1 



equations with 'energy dependent potentials' due to the analogy with the 
Schrodinger equation in Quantum Mechanical context, whose spectrum rep- 
resents the energy levels of the Quantum Mechanical system. Some of these 
integrable systems appear as water waves models, most notably the Kaup- 
Boussinesq equation [221 EOJ [7] and the two-component Camassa-Holm equa- 
tion [ini E21 UM ■ Other systems of this type are studied e.g. in El US] • 
In what follows we study an integrable system which arises as a compat- 
ibility condition of the following two linear operators (Lax pair): 



* XX 



-A 2 + Xu(x, t) + -u 2 (x, t) + T](x, t)J ^0, t] A) (1) 
*t = - (a + ^u(x, t)j if> x (x, t- A) + ^u x {x, t)^(x, t- A). (2) 

Here k is an arbitrary constant, while A is the spectral parameter. The con- 
sistency condition ip xx t(x, t; A) = iptxx(%, t\ A) produces a system of equations 
for the functions u(x, t) and r)(x,t): 



u t + r] x + y- + kJ uu x = 0, (3) 

Vt ~ ^u xxx + (urj) x - Q + kJ ur] x - k ^ + u 2 u x = 0. (4) 

Upon choosing k — — 5 these simplify to the well known Kaup-Boussinesq 
(or KB for short) equation: 

Ut+T)x + uu x = 0, (5) 

Vt ~ \u xxx + (ur]) x = 0. (6) 



The KB equation is introduced as a water-wave model in [22] where also 
the inverse scattering is studied for functions with constant limits at x — > 
±00. As a water-wave model it also appears in [301 HHJ E EH], the hierarchy 
of Hamiltonian structurs is given in [27], specific solutions are studied in 
[21~1 El El]. Energy- dependant spectral problems like ([1]) are studied also in 

[nanu [m[i9i[2a m\ m\. 

Our aim will be to formulate the inverse scattering as a Riemann-Hilbert 
Problem (RHP) in the case when u(x, t) and rj(x, t) are real, rapidly decaying 
functions at x — > ±00, taking into account the underlying reductions and to 
obtain the simplest soliton solutions. 



2 



2 The spectral problem 

Introducing an auxiliary function 

w(x,t) = -u 2 (x,t) + 7](x,t) (7) 

we consider the following two 'conjugate' spectral problems related to ([1]): 

V xx (x, A, cr) = (-A 2 + a\u(x) + w(x)) $>(x, A, a), (8) 

where a — ±1. The t-dependence will be suppressed where possible for the 
sake of simplicity. 

We specify that u(x),w(x) as well as rj(x) belong the the Schwartz class 
of functions (the space of rapidly decreasing functions) iS(lR). It follows from 
this requirement, that solutions i[)i(x,t; A) and i[) 2 (x,t; A) exist such that, 



ipi (x, X, a) — > e iXx 
ip 2 (x, A, cr) ->■ e +iXx 



x -> +oo (9) 



Similarly we define a basis of eigenfunctions for (JSJ) according to 

x ->■ -oo (10) 



0i (x, A, cr) -> e iAx 
c/> 2 (x,A,cr) e +iAl 



These eigenfunctions are called Jost Solutions. Since the Jost solutions 
oscillate when A 6 R, the real spectrum fills in the real line. 

The bases {ipi(x, A, cr), ip 2 (x, A, cr)} and {c/>i(x, A, cr), c/> 2 (x, A, cr)} consti- 
tute independent bases of solutions to (|S]) and as such, we may write 

(pi(x,\,a) \ = f T n (\,a) T 12 (X,a) \ / ifa(x, A, a) , 
2 (^,A,cr) y V T 2 i(A,cr) T 22 (A,(t) y \ V 2 (x, A, a) 

The matrix 

Tn(A,(x) Ti 2 (A,a 



T ^> = U»W £m) (12) 

is the scattering matrix for spectral problem (jSJ). 

Under the involution (A, a) — >■ (— A,— cr), the potential in (jHJ) remains 
invariant. Therefore the eigenfunctions ^(^ A, cr) and ip(x, —A, — cr) are solu- 
tions to the same spectral problem. Since the asymptotics of these solutions 
do not depend on cr, it follows that 

ipi(x, A,cr) = i/) 2 (x,-\-<r) 

(f>x(x,X,a) = <f) 2 (x, -A, -cr), (13) 
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Thus, we can write the two bases using just one of the functions, say ip(x, A, a) = 
ipi(x,X,a) and <p(x,X,a) = <f)x(x,\,a) as ip2(x, A, cr) = ip(x, — A, — a) and 
02 A, cr) = 4>{x, -A, -cr). 

When u(x,t) and r](x,t) are real, the spectral problem (jHJ) is invariant 
under Z 2 reduction group [25], i.e. it has the following property: if if)(x, A, a) 
is an eigenfunction, so is ^(x,X, a). Comparing the asymptotics again, we 
conclude that this coinsides with the second Jost solution, i.e. 

$(x,\,a) = 'tl>(x,-\,-(T) (14) 

Thus, for A 6 R we also have ip(x, A, a) = ip(x, —A, —cr). From this, and 
ffTl~j) it follows that the scattering matrix T(A) may be written in the form 

(i5) 

for spectral parameter A G R. 

We now have the following relationship between 0(:r, A, a) and the Jost 
solutions i/j(x, A, cr), ^(x, A, a), 

0(x, A, cr) = a(A, a)if)(x, A, cr) + 6(A, a)ip(x, —A, — cr). (16) 

Furthermore, for any pair of solutions f\ and / 2 to ([HD the Wronskian of the 
pair is independent of x, 

d x W[h, / 2 ] = ^(/x^/2 - WO = 0. 

In particular, it follows that the Jost solutions satisfy the following condition 

W[(f){x, A, a), <j>(x, A, a)] = W{ip(x, A, a), ^(x, A, a)} = 2i\, (17) 

which clearly follows from the asymptotic behaviour of {ip(x, A, a)ip(x, A, cr)} 
and {0(x, A, cr), 0(x, A, cr)} as \x\ — > oo. It follows from ffT6l) and ffTTl) that 

detT(A,cr) = |a(A,cr)| 2 - |6(A,cr)| 2 = 1, A G E. (18) 

3 Asymptotic behaviour of the Jost solutions 

Since the functions u(x) and w(x) are Schwartz class it follows that the 
solution ip(x, A, cr) have asymptotic behaviour such that 

if) xx {x, A, a) ->■ -\ 2 e~ iXx , x -)> +oo. (19) 
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Consequently, we make the following ansatz for the asymptotic expansion as 
|A| — > oo, 



if>(x,X,a) 



X (x,a) + jX 1 (x,a) + O(\- 2 ) 



-i\x 



(20) 



where the function Xq(x,<j) and X\(x,a) behave asymptotically according 
to 



X (x, a) — > 1 
X 1 (x, cr) 

The substitution of fl20|) into ([8]) gives 

^X (x,cr) i 
Xo(x s< 7) ~ 2 
au(x)Xi (x, cr) + 2i<9 x X 1 (x, cr) 



x — >■ +oo. 



au(x) 



(21) 



(22) 



-X (x,cr) + d 2 x X (x,a). (23) 



Using the conditions in (j2ll) we may easily solve (122ft . (|23i) to give the fol- 
lowing expressions for Xq(x, a) and Xl(x, cr): 



X (x,Cr) 



exp 



Xi(x,a) = X Q (x,a) 
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u(x')dx' 



-u[x\ 



u 2 (x') + 4w{£))dx' 



Similarly, we obtain analogous expressions for ip(x, A, cr), i.e. 



e -»(Aar+f f™u(x')dx') 



ip(x,X,a) 
<P(x,X,(t) 

where the functions £i(x) and (i{x) are given by 



1 + -£ 1 (x,<t)) + ... 



(24) 

(25) 
(26) 



£i{x,a) 
(i(x,a) 



a 



u(x) 



ju(x) 



u (x') + 4w(x'))dx' , 
u 2 (x') + Aw(x'))dx'. 



(27) 



4 Analytic behaviour of the Jost solutions 



We now define the related function 

e 



(x, A, a) = e lXx (/)(x, X, a) ->■ 1, x -)■ -oo. 



(28) 



Using 

e iXx cf> x {x, A, a) = xi +) (x, A, a) - iX X (+) (x, A, a) 
along with the spectral problem in (jHJ), we may write 

X { x + J(x, A, a) = {Xau{x) + w(x)) X i+ \x, A, a) + 2i\ X { x +) (x, A, a). (29) 



Meanwhile the asymptotic expansion in A appearing in f l26|) suggests the 
following integral representation for x^ + K x , A, a), 

/x 2i\(x~x') _ i 
— P(x', A, a) X {+) (x', A, (30) 
2tX 

for some P(x, A, a) G iS(lR). From this integral representation differentiating 
twice we obtain 

X&\x,\,(t) = P(x,X,a)x i+) (x,X,a) + 2tXxi +) (x,X,a). 

From (I29p and (131 p we determine 

P(x, X, cr) = Xau(x) + w(x), (31) 

and so we may write 

-^—(Xau(x f ) + w(x'))x {+ \x f ,X,a)dx>. (32) 

Of particular importance and clear from (13"2"|) is the analytic properties of 
X^(a;, A, a). We can see that for all values of x the kernel of the integral 
above is finite for all values of A such that Im A > 0. Therefore X ( x i \ a ) 
and <j)(x, A, a) are analytic in the upper half plane C+. It obviously follows 
that x ( x 7 K a ) i s analytic for A G C_. 
In a similar manner we may define 

X { ~\x,X,a) = e iXx ip(x,X,a) ^ I, x -> +oo (33) 

from which it follows that 



oo e . 



X { -\x,X,a) = l- I e2lX(X 2 ^ - (Xau(x)+w{x))x H (x,X,<7)- (34) 



It is immediately clear from ( 15%|) that ^ ^(x,A,cr) and therefore ip(x,X,cr) 
are analytic throughout C_. 
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Next we introduce new notation for later convenience, 

rx /"OO 

U-(x) — - I u(x')dx' and 0J+(x) — — / u(x')dx' . (35) 
With this we may rewrite (125|) . (1261) as follows 

A, a) = A, ff ) e <C>-W»)) = 1 + I^x, a); 

A, a) = <p{x, A, (**-«"-(*)) = i + l&fa a ). ( 36 ) 

A 

To obtain the analytic properties of ip(x, A, a) and 0(x, A, a), we note that 

$(x, A, a) = x ( ~ } (^ A, t^e*"*^ and 0(x, A, a) = x {+) (^ \ a)e~ iff ^ (x) . 

Since u(x) is Schwartz class and independent of A and given the analyticity 
of x ( x > K a ) throughout C± respectively, it follows that 0(x, A, a) and 
ip(x, A,cr) are also analytic throughout C + and C_ respectively. 

5 The ^-dependence of the scattering data 

We may rewrite the second member of the Lax pair in terms of the auxiliary 
function u(x) and add an arbitrary constant 7, without effecting the physical 
equations of motion, to obtain, 

V t {x, A, a) = - (a\ + ^u(x)^j ty x (x, A, a) + ^7 + ^u x {x)J A, a). (37) 

In particular we may write 

<f) t (x, A, a) = - ^crA + ^m(x)^ (f) x (x, A, a) + ^7 + ^(2;) J 0(x, A, cr). (38) 

However, we also note that along the discrete spectrum we have the scat- 
tering relation ( TT6|) . from which we may obtain the asymptotic behavior of 
4>t(x, A, cr) as x — > +00, namely 

<j> t {x, A, a) ->• a t (X, a)e~ iXx + b t (X, a)e +lXx . (39) 

Using the r.h.s of ( 138]) along with ( !T6|) . we find as x — > +00 that 

0i(x,A,cr)^ - (TA[-iAa(A,cr)e- iAx + iA6(A,cr)e iAx ] 

+ -flaiX,^-^ + b{\,a)e +iXx ], (40) 
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where we have made use of the fact that u(x) is Schwartz class and vanishes 
when x — > ±00. Making the choice 7 = — iaX 2 , the t-derivative of a(A, a) 
vanishes. It follows that we may write 



a t {\,cr) = a(X,a,t) = a(X,a,0), , 
6 t (A, a) = -2tX 2 b(X, a) b(X, a, t) = b(X, a, 0)e~ 2i<TA2 * ' 



Along the discrete spectrum, we have a(X n ,a) = and therefore instead 
of fl39l) we have 

^(x,A n ,a)^6„,i((T)e iA ^. (42) 

Instead of (14"U|) we have 

Mx, X n , a) ^ -aX n [iX n b n (a)e iX ^} + 7 (A n )6 n (a)e iA " a; , (43) 

and thus 

b n!t {a) = -2iaX 2 n b n {a), (44) 

from where 

b n (a,t) = b n (a,0)e- 2t ^ t . (45) 

6 Conservation Laws 

We may derive a collection of conserved quantities from the spectral problem 
introduced in the previous section in ([8]). To proceed we first introduce the 
function 

<£? x (x,X,(j) 

Pt?M= Hx ^ a y (46) 
Differentiating once with respect to x we find 

p 2 + Px = -X 2 + aXu + w. (47) 

Using this result, along with the Lax pair in ([8]) and (1371) we find upon 
differentiating (T4*6l) with respect to t that 

Pt = \ («x - 2pw - AaXp) x . (48) 
Using the the fact u(x) and w(x) are Schwartz class, we see from fT48|) that 



+00 



p t (x, X,a)dx = 0, (49) 
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that is to say 

J (A) = / p(x,X,a)dx (50) 



is a generating function for the conserved quantities. We may expand it in a 
power series in A according to 



X = XL, +X 2 + =| + =| + 0(\- 3 ), ill ) 



where Xx,X 2 , etc is an infinite sequence of conserved quantities. Next, we 
expand p(x, t, A) as a power series in A 



p(x, A, a) = a + p (x) + + ^ + 0(X-% (52) 

and use it in (T4"Tj) . then the terms of equivalent order in A give 



PoW = (53) 



and result of (ED) it follows that 



«7 / " +0C 



OO t/ — oo 



2^2 = / po(x)dx = —— J u(x)dx. (54) 



So we see that 



1 



ai = - u(x)dx (55) 

^ i-oo 

is an integral of motion. Following a similar procedure, we find the next 
conserved quantities to be 

hoo 



J 3 = —= / (u 2 (x) +4w(a;))ofe, 



X 4 = / n(x)(M 2 (x) + 4w(x))(ia;. (56) 

!6 J-oc 

One may continue a process of iteration indefinitely, whereby an infinite 
series of such conserved quantities is generated from the u(x) and w(x), and 
therefore from the physical variables u(x) and r)(x). 
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7 Analytic continuation of a(A, a) 



Returning to ([16]) we see that we may re-write the scattering coefficient 
a(X, a) in terms of the ^-independent Wronsikian, 

iy[0(x,A,a),^(a:,-A,-a)] 
a(A, a) — . (57) 

Since the two eigenfunctions in (l5"Tj) are analytic for A G C+, a(A,cr) allows 
an analytic continuation in the upper half complex plane. From (15?]) with 
(I25p - 0261) we obtain the asymptotic behavior of the scattering coefficient, 

lim a(A, a) = e iaai (58) 

A|— >oo 

where ot\ is the conserved quantity (155}) . We make further the assumption 
that a(A, a) has a finite number of simple zeros A n G C + , n = 1, 2, 3, . . . , N. 
We introduce the auxiliary function 

N A- A 

A(A,a) = e— i nW ia ( A ' cr )' ( 59 ) 

which is analytic without zeroes in C+. It follows from (159]) that 

\A(\,a)\ = \a(\,a)\i A G R, (60) 
Next, we also see from (1581) and (1591) that 



lim lnA(A,o-) = 0, (61) 

A|— >oo 

and so, ln^4(A,a) is analytic throughout C + and vanishes as |A| — > oo. 
We also have from (1601) 



In A(X, a) = In \A(X, cr) \ + i arg A(A, a) = In |a(A, cr) | + i arg A(A, a) 
for A G M. We make use of the Kramers- Kronig dispersion relations, 

oo 

In |a(A,cr)| = - f dX 



71 X'-X 



-oo 
oo 



1 f l n l a (A', cr) I „ . . . 

argA(A,a) = -- j \,_\ n dX', (62) 
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for A G R, where the dashed integral denoted the principal value part of the 
integral. Then with (162 p we have 



OG 



\nA(X,a) = ln\a(\,v)\-tJ hMXlpl dX > 

7T ./ A' — A 



oc 



I f ln|a(A'^)| dA , AeR _ (63) 



-oo 



A — A r 



m .1 A' - A - i0+ 

Meanwhile, f loTJj) gives 

In A(A, o-) = -taai - V * - ^ + In a(A, a). (64) 

n=i X ~ Xn 

Using (16 3 p and flMj) . we find that for real values of A we may write 

N 00 

x • A — A n 1 /" lna(A,c)l „. . . 

In a (X.a )= %oa x + > — + — / ' dA ; 65 

v ' ' ^ A - A n vri 7 A' - A - i0+ v ; 

— oo 

and for A G C + the analytical continuation is 



A' 



A — A n 1 f In |a(A', cr)| , 



In a(A, a) = + V + - / '"^ \ U)l dX>. (66) 

/ — ' A — \ n Til J A — A 



n=l 



-oo 



8 The Riemann-Hilbert problem 



We may re-write the expression fTl6|) in terms of the new analytic functions 
<f)(x, A, a),ip(x, A, a) using (fT4"|) for A G R, as follows 

<^>(x A crW " 01 

' ' = ^{x,Xa)+r{X,a)^{x,X,a)e 2i{Xx+au,+ ( - x) \ (67) 



a(A, cr) 

r(A,cr) = 6(A, a)/a(A, cr). The function — 'A s — is analytic for Im A > 0, 
while ^(rr, Act) is analytic for Im A < 0. Thus, equation fl67|) represents an 
additive Riemann-Hilbert Problem (RHP) with a jump on the real line, given 
by 

r(A,a)^(x,A,a)e 2i(Ax+CTaJ+(x)) 
and a normalization condition limiA^oo ip(x, A, a) = X (x, a). 
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Figure 1: The integration contours C + and C~ are the closed paths in the 
upper and lower half planes correspondingly; T± are the semicircles with an 
infinite radius. 



In this section we will follow the standard technique for solving RHP. 
We integrate the two analytic functions with respect to ^ ttzt( - ) over the 
boundary of their analyticity domains, using the normalization condition. 
In our case the domains (the upper C + and the lower C_ complex half- 
planes) have the real line as a common boundary and there we relate the 
integrals using the jump condition. The RHP approach for various equation 
is presented in PH [231 E21 EH] • 

We now choose some A G C_ and integrate the left-hand side as follows, 

1 f 0(x,\',a)e iaai , * ^ n) (x,a)e iaai 
2ri f c+ a(X',a)-(X'-X) dX = ^ a n (a) ■ (A n - A) (68) 

where C + is the contour in the upper half plane shown in Fig. 1, 
a n (a) = (j^^j ^ 0, tf n \x, a) ee A n , a). 

We may write the integral as such because A e C_ and so -j^y is ana- 
lytic throughout C + . Furthermore, a(A) is analytic with finite number of 
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simple zeros, A„ in C+, and the function (f)(x, A) is analytic throughout C H 
Alternatively we may expand the integral as follows 

1 £ 0(x,A»e iQl 

~~ -aX = 



2m f c + a(X',a) ■ (A' - A) 

1 f°° 6(x, A', a)e iai 1 f (f)(x,X',a)e iai 

~ -dX' + — / =± \ ' rrdX'. (69) 



2m y_ 00 a(A') ■ (A' - A) 2m J T+ a(A', a) ■ (A' - A) 

Using the asymptotic properties of a(A,cx) and A, a) along with the re- 
lationship f )67|) . we find 

" (w) (x,a)e^ J_ r 1 4 _L r #^ A >) ,y 

^ d(A„) • (A n - A) 2mJ r+ X'-X + 2mJ_ 00 A' - A 

1 Z" 00 r(AV)e 2i ( A ' x+ffW +^V^, A',cr) 

+ / ; ; — = dX . 

2m A' - A 

(70) 

Next we obtain an expression for the line-integral 

1 f+°°±(x,X',a) 

-dX , 



2m J_ 00 A' — A 

by considering the integral over the contour C~ , shown in Fig. 1. Since 
A G C_ and ip(x, A, a) is analytic therein. In addition the contour C~ is 
clockwise, so it follows that 

1 / ib(x,X,a) , 

sslfe^" (71) 

Expanding the integral, we have 

1 f i/;(x, X',a) 1 f +oc ip(x, A', a) 1 f ^(x,X',a) 

2mf c -^y^r dX= 2mJ 0O =X^X~ dX + 2m J T = A 7 ^A _ ' 

( (72) 

Using the asymptotic properties of i/j(x, X,<t) as |A| — > oo with A G C , we 
have, 



l r°° ib(x,X'(j) i r i 

ii =1J — ——dX = -ip(x, A, a) — / -^—dX' (73) 

27ri 7.^ A' - A - v ' ' ; 2m J r A' — A v ; 



We can also make use of the following result when it comes to substituting 
this expression in (170|) . 



' dX' - [ —^—dX' = 2m. (74) 



r+ — X Jr— X' — A 
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Upon making these substitutions we find the following integral representation 
for ij)(x, X, cr), A G C_: 

(j) (n) (x,a)e iai 1 f°° r(X')e 2l ^' x+aul + {x ^i)(x,X',a) 

ii>(x, A, cr) = 1— > = . . . , r-H / : r - = dX . 

~ ' ^d n (a)(X n -X) 2m X'-X 

(75) 

Since at the points of the discrete spectrum <fi(x, X n , cr) = b n (a)ip(x, X n , a) we 
have 

^ n) (x,a)e tai 



a n {a) 

where we define 

RJ<J 



iR n (a)e 2l ^ M + (x »l(x, X n , a) (76) 



ia n {a) 



The Riemann-Hilbert problem is reduced to the linear singular integral equa- 
tion for ip(x, A, a) 



N 



±{x,X,a) = l-i^2 



n=X 



R n (a)ij_(x,X n ,a) 
(A~^A) 



2i(\ n x+ou> + (x)) 



1 r°° r(X',a)ib(x, X'.a) , , „ , 
2m J ^ X'-X 1 ; 

In addition to (1771) we have an analogous system written at the points 
A = X p G C_, p= 1,2,..., N: 



E N R n (a)ib(x, X n ,a) „.,, , ^ 

(A n — A p ) 



n=l 



2~ki ,/_ 00 A' — A p 

Finally, the fact that at A = the Jost solution i/j(x, 0, cr) does not depend 
on cr gives ip(x, 0, a) = if)(x, 0, — cr) or an algebraic system for e 2i ™+(x) : 

n- ci ^(x,0, cr) ib(x,0,a) 
e 2iou> + (x) _ za ; ; ^ __ z.v ' ; / c^ 

^(x,0,-a) ^(z, 0, cr) ' 

The system (|77jl . fl78|) . ( 179]) allows for the determination of both the Jost 
solution and the potential functions of the spectral problem in terms of the 
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scattering data. Note that the time- dependence of the scattering data is 
known from flUD, 

r(X,a,t) = r(\,a,0)e- 2l ° x2t , R n (a,t) = R n (a,0)e- 21 ^ 1 . (80) 
Thus, the complete set of scattering data is 

r(A,M), A„, R n (a,0) (n = 1, 2, . . . , N). (81) 

Also, it is sufficient to know the scattering data for a = 1, because of the 
%2 involution, which holds on the scattering data too: 

r(A, -a) = f(-A, a), R n (-a) = R n (a). (82) 

9 Reflectionless potentials and soliton solu- 
tions 

The so-called reflectionless potentials are a subclass which corresponds to a 
restricted set of scattering data: r(A, a) = 0; A G R. Then the system (f77|) . 
(ITS"]) . ( ITUj) is algebraic, and the solutions of the PDE are called solitons. 

The simplest case is the N = 1-soliton solution, so we start first with this 
case. From ( 177j) we have 

Ri (a)%b(x, Ai, a) „.,. , 
^(x; A, cr) = 1 - z V [- L e M*i*+™+i*)) } \ e c_. (83) 

Ai — A 

We notice that ^(x, A, a) has an unique pole at Ai and ip(x, —A, — cr) has 
an unique pole at —Ax. Due to ( flij) these two poles coincide, i.e. Ai = — Aj 
and therefore Ai = zV is purely imaginary, v > is real. 

Solving for ?/>(x, Ai, cr) we find 

^ ") = ; rfeMig^ — (84) 

1 + 4Af 

Then ( 18 3 p takes the form 

^(x,A,a) = l + T ^ ^ ^ (85) 
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Furthermore, we can relate the real and imaginary parts of the complex 
constant to two new constants, say Xq and t as follows: 



2u 



g4rao — 2iaf 2 to 



Now ip(x,\,a) in fl85l) depends only on x — x , t — t and due to the 
translational invariance of the problem, without loss of generality, we can 
choose x = and t = 0. This simplifies f lSo^) to 



^{x, A, a) = 1 + ^— (86) 



2,ilS ^,—2h'x+2iau t+2iaui+(x) g— Avx 

A — iv 1 — e~ 

Then ([79]) gives 

i _j_ 9 2vx— 2iav 2 t _j_ _— 4^3; 
J _|_ 2g — 2ra+2«cw 2 t _|_ g— 4ra 

From a; + (x,t) and fl35]) we can recover u(x,t): 

sin(2z/ 2 t) sinh(2z/a;) 



t) = v- 



cosh 4 (i/a;) cos 2 (z/ 2 t) + sinh 4 (z/x) sin 2 (^ 2 t) 
On the other hand, we also have 



Tp(x,\,a) = 1 + - 



4 

which can be compared to (186]) : 



OG 



a n(x) + I I {u 2 {x) + 4w(a/))da/ 



O 



Oin p—2vx+2iau 2 t+2iaLOJ r {x,t) „— 4vx / i 

*(*A») = l + x- ! i- e -~ + ° (^)- <*" 



Since u + (x,t) and u(x,t) are already known, we find w(x,t) and finally 
T](x,t). With (]7j) we compute 

M 2 + 4w = 2(k + ^)m 2 + 4r?. 

For the Kaup-Boussinesq and 

+ 4w = 4ri = -Ad 2 x In [(1 + e" 2ra ) 4 + (1 - e~ 2ra ) 4 tan 2 uH] , 



u 2 
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Figure 2: Snapshots of the solutions of the KB equation (IHHj) . (190]) for three 
values of t. The first panel is before, the third panel is after the blowup. 



2 cosh 6 (z/x) cos 4 (z/ 2 t) + f sin 2 (2z/ 2 t) sinh 2 (2z/x) -sinh 6 (z/x) sin 4 (z/ 2 t) 

^ — — 2 z/ , 

[cosh 4 (z/x) cos 2 (z/ 2 t)+sinh 4 (z/x) sin 2 (z/ 2 £)] 2 

(90) 

The solution ( |88|) . (190]) is presented on Fig. |2j Note that w is an odd and 
i] is an even function of x. The solution is of 'breather' type and develops 
singularities 'infinitely' close to x = at countably many isolated values of 
t. 

The next case is a solution with N = 2 discrete eigenvalues. Due to ( II 4ft 
there are the following situations: 

(i) Both eigenvalues are on the imaginary axis: Ai = iu%, A2 = iv-i for 
some real and positive v\ and z/ 2 ; 

(ii) A 2 = — Ai, R%{p) = Ri(—a). For the (ii) case from ( 177|) we have 



ip(x; A, a) 



1+ie 



i2i(ff)e 2iAl!B ^(a:, Ai, a) R^-v^^^ix, Ai,^r) 



A-Ai 



A + Ai 



(91) 
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From (19Tj) we obtain a linear system of four equations for the quantities 
ip(x, Ai, ±cr) and their complex conjugates by writing (l§Tj) for A = Ai, the 
same with a replaced by —a and their complex conjugates. 

The case with N > 2 eigenvalues is always a combination between (i) 
and (ii) - in general it involves eigenvalues on the imaginary avis as well as 
conjugate couples A& and — A&. 

10 Conclusions 

We have outlined the inverse scattering for the spectral problems of the form 
(jBD with real functions in the potential, which necesitates the Z 2 reduction 
(I14p . The soliton solution in the case of a single pole of the eigenfunction does 
not have the form of a travelling wave and develops singularities with time. 
This solution is probably not relevant for the theory of water waves. There 
is another feature of this type of equations which points in the direction that 
the purely soliton solutions are probably not the ones which are observed 
in the context of water waves. Indeed, since rj is the deviation from the 
equilibrium surface, then one expects that its space-average value is zero, 
f_ r](x,t)dx = 0. However, the trace identities which can be derived easily 
(see e.g. [23]) for the iV-soliton solution of the KB equation lead to the 
following result: 

/OO /"OO N 

r](x,t)dx = - / (m 2 + 4w)dx = — 4 ImAfc. 
-oo 4 J —oo k=1 

By assumption ImAfc > since A& are in the upper half complex plane. 
Thus, we have the following 'mostly negative' result for the iV-soliton solu- 
tion: 

/oo 
r](x, t)dx < 0. 
-oo 

This results indicates that the water wave solutions are related only to the 
continuous spectrum and are therefore unstable. This agrees with the fact 
that the travelling wave solutions to the Euler's equation with zero surface 
tension are unstable. 
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